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Abstract. We give a construction of a self-similar tiling of the 
plane with any prescribed expansion coefficient A G C (satisfying 
the necessary algebraic condition of being a complex Perron num- 
ber). 

For any integer in > 1 we show that there exists a self-similar 
tiling with 27r/m-rotational symmetry group and expansion A if 
and only if either A or Ae 27 ™/" 1 is a complex Perron number for 
which e 2 ™/" 1 is in Q[A], respectively Q[Xe 2m / m }. 



1. Introduction 

By a tile we will mean a compact subset of R 2 which is the closure 
of its interior. A tiling of the plane is a collection of tiles T = {Tj}j £ /, 
whose union is R 2 and which have pairwise disjoint interiors. 



In ||18|| , Thurston introduced the notion of self-similar tiling. This 
definition was motivated from several sources: most notably, in tilings 
arising from Markov partitions for hyperbolic toral automorphisms, 
one-dimensional tilings arising from substitutions, and also the Penrose 
tilings, which have a simple subdivision rule to create new tilings out 
of old. 

A self-similar tiling is a tiling T = {T{\i & i of the plane which has 
the following properties: 

1. There is an equivalence relation ~ on tiles, with a finite number 
of equivalence classes, such that Tj ~ Tj implies Tj is a translate 
of T t . 

2. There is a homothety ip: C — > C, <p(z) = \z, such that the image 
of a tile is equal to a union of tiles of T, 

3. If Tj ~ Tj, and Tj + c = Tj, then for each tile Tk C f{Ti), y?(c) +Tfc 
is a tile of T equivalent to Tk- 

4. The tiling is quasiperiodic (see definition below). 

An example is shown in Figure [T] (more information about this example 
is given in section H). 

Let us make a few remarks regarding this definition. 
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Figure 1. Part of a self-similar tiling (see section 

1. The third condition simply states that two equivalent tiles "subdi- 
vide" in the same way. 

2. Here by quasiperiodic we mean, any "arrangement" (finite col- 
lection of tiles, up to translation) of tiles in the tiling reoccurs at a 
bounded distance from any point in the tiling. More formally, for 
any r > there is an R > such that for all x, y E C a translate 
of the arrangement of tiles in B r (x) occurs as a subset of in Bn(y). 
This property is almost the same as the local isomorphism property of 
Radin and Wolff HT51, and is identical if one assumes in their case that 



there are a finite number of local arrangements (which in our case is a 
consequence of the definition). This condition is really a compactness 
condition: indeed, this condition is equivalent to the condition that the 
orbit closure Y of the tiling (under the action of R 2 by translations) 
is minimal: every R 2 -leaf of Y is dense in Y [15]. Under reasonable 



assumptions on the subdivision rule for tiles, this condition is easily 
checkable, see Lemma |], below. 

3. By replacing "translations" in the first and third conditions by some 
other group of homeomorphisms of M 2 , e.g. by isometries, one arrives 
at a more general definition of self-similar tiling. So technically we 
should refer to our definition as that of a translation-self-simil&r tiling. 
However for succinctness (and to agree with Thurston's terminology) in 



THE CONSTRUCTION OF SELF-SIMILAR TILINGS 



3 



this paper we shall simply say "self-similar tiling" , with the translations 
understood. The results as stated below are unknown for isometry-self- 
similar tilings; for similarity-self-similar tilings, see |ll] . 



Thurston proved that the expansion coefficient A of a self-similar 
tiling had to be a special type of algebraic integer, a complex Perron 
number. He indicated, without proof, that for each complex Perron 
number there existed a tiling. In this article we give a construction of 
a self-similar tiling for any complex Perron number. 

A complex Perron number is an algebraic integer which is strictly 
larger in modulus than its Galois conjugates (the other roots of its 
minimal polynomial), except for A, its complex conjugate. 

Theorem 1 (Thurston flT8|). The expansion of a self-similar tiling is 
a complex Perron number. 

In section |] we will prove by construction: 

Theorem 2. For each complex Perron number A there is a self-similar 
tiling with expansion A. 

The methods we use are very geometric and are inspired by conver- 
sations with and ideas of Thurston. 

This construction allows an easy corollary, whose proof we give in 
section ^ 

Theorem 3. A is the expansion of a self-similar tiling invariant under 
a rotation e 27 ™/ n (n an integer) if and only if either A ^ K. and A is a 
complex Perron number for which e 27 ™/ n £ Q[A], or \e 2m / n is a complex 
Perron number and e 27 ™/ n E Q[Ae 27n//n l. 



Note that in the case n = 2 we conclude: any complex Perron number 
is the expansion for a self-similar tiling invariant under x — > —x. 

For an example of this theorem, consider the symmetric Penrose 
tiling, which is invariant under rotation by 2tt/5. One possible expan- 
sion of this tiling is the degree-4 complex Perron number A = 1 + e 27 ™/ 5 
and Q[A] = Q[e™/ 5 ]. 

In section |6] we will give a simpler construction for self-similar tilings 
whose expansion coefficients are of a particular type: those non-real 
complex Perron numbers A satisfying an equation of the form 

A n - aX n - 1 + b\ + c, 

where a, b > 0, c > 1 are integers. This construction has the advantage 
of being very concrete. 
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2. Historical remarks 

Self-similar tilings appear in several quite different contexts. In fact 
up to now they have not really had a life of their own, but have rather 
been studied independently by people working in different fields. 

Self-similar tilings did not really exist before the computer age, al- 
though some precursors can be found 0. One could also trace back 
the roots to the work of Adler and Weiss ][]], Bowen || and Sinai |16| 



on the construction of Markov partitions, which can be interpreted as 
the first general constructions for self-similar tilings. 

Most of the first explicit examples of self-similar tilings were curiosi- 
ties, found in popular books on fractals with fanciful names such as the 
"Heighway dragon curve" and "Gosper flowsnake" . 

The Penrose tilings were another curiosity until some nonperi- 
odic crystals displaying 5-fold "symmetry" were actually discovered in 
nature; these "quasicrystals" gave life to a new field of crystallography. 

The theory of wavelets is another new field which has caused peo- 
ple to become interested in "self-replicating" tilings (self-similar tilings 
with one tile type), and much work has been done since on this special 
class of self-similar tilings 0, |9|, |17| ^ . 



Our current construction has its origins in another point of view. 
Working in the field of symbolic dynamics, Lind |13| gave a characteri- 



zation of the largest eigenvalues of nonnegative integer matrices, which 
gives as corollary a characterization of the expansions of self-similar 
tilings of the line. 



Then Thurston |TJ| in 1989 took a step up from this result, defining 
self-similar tilings and proving the neat Theorem [I] above. As a thesis 
student of Thurston, I extended his result to n-dimensional self-affine 
tilings (in which the expansion is a linear map of M. n diagonalizable 
over C) Jig. 

Not much later Praggastis gave a more explicit construction of 
Markov partitions for (certain) hyperbolic toral automorphisms using 
the theory of self-similar tilings. This work is also very close to that 



of Vershik [19[ who constructs such Markov partitions arithmetically, 
that is, using radix representations. 

Self-similar subdivisions (of various sorts) occur nowadays in the the- 
ory of rational maps, Kleinian and hyperbolic groups (via the action on 
the sphere at infinity or space of geodesies using an automatic struc- 
ture), wavelets, crystallography (in quasicrystals), symbolic dynamics 
and number theory (in radix representations), and logic (in Wang tiles, 
aperiodic tilings, and computability: surprisingly, the only known ape- 
riodic tilings are built from self-similar patterns). 
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3. Background 

3.1. the subdivision rule. The tile types in a self-similar tiling are 
the distinct equivalence classes of tiles. To a self-similar tiling with n 
tile types is associated a nonnegative n x n integer matrix M = (rriij), 
the subdivision matrix, in which rriij is the number of tiles of type 
j that a tile of type i subdivides into upon scaling by A. By property 
3 this matrix is well-defined. 

The Perron (i.e. largest) eigenvalue of this matrix is AA, whose cor- 
responding eigenvector is the vector of tile areas. 

Recall that a nonnegative matrix is called primitive if some power 
is strictly positive. In the case of a self-similar tiling, the subdivision 
matrix must be primitive, since every tile, when scaled by A n for large 
enough n, must contain a copy of every tile type by quasiperiodic- 
ity Primitivity also implies that the spectral radius AA is a Perron 
number, that is, AA is a real algebraic integer strictly larger than the 
modulus of any Galois conjugate. 

We have the following result (whose proof is an application of the 
definitions) which we shall use later. 

Lemma 4 ([14|], see also JTTJ ) . A tiling which satisfies the first three 
hypotheses of a self- similar tiling, has a tile with the origin in its inte- 
rior and has primitive subdivision matrix, is quasiperiodic, and hence 
self- similar. 

3.2. The Delauney triangulation. We review here the definition of 
the Delauney triangulation. Given a discrete set of points P in M 2 , 
the Delauney triangulation is the triangulation of the convex hull of 
P, with triangles having vertices in P, with the property that, for 
any triangle t = (t>i, t>2, fs), the circumcircle of t contains no point 
of P. Such a triangulation exists (for example, one can see this by 
using stereographic projection of the plane to the sphere in M 3 ; the 
convex hull in M 3 of the image of P is a polyhedron whose faces are 
(generically) the triangles in the Delauney triangulation) and is unique 
unless some four or more points lie on the same empty circle. 

Another property that we will use is the local nature of the triangu- 
lation: if for some small e there is a point of P within a distance e of 
any point of P>i(x), then the set of triangles with vertex x depend only 
on those points in B 3t (x). (To see this, notice that an empty circle 
passing through x has diameter at most 2e.) 



4. Proof of Theorem |2] 
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4.1. sketch. The idea of the construction is as follows. We define the 
set T of "archtypes" of the tiles in our eventual tiling to be the set of all 
triangles with vertices in an appropriate lattice and edges of bounded 
length. When we multiply a triangle in T by X n for large enough n, 
we can subdivide it approximately into triangles in T. (This is where 
we use the complex-Perron-number property of A). However we need 
a way to define the subdivision of triangles so that, if two triangles 
are adjacent, their subdivisions agree at their common boundary. So 
rather than just keep track of a triangle we keep track of a triangle and 
all its immediate neighbors. This gives a new larger set of archtypes 
X. The subdivisions of elements of X can now be defined so that they 
agree whenever there is an overlap between them. 

One small hitch is to make this subdivision rule primitive. To ensure 
this we define a special tile T that occurs in the subdivision of any tile, 
including itself. Our set of tiles will then be the strongly connected 
component of the "subdivision graph" which contains To. 

Finally we redraw the boundaries of the elements of X U {To} by a 
recursive process so that the subdivision is exact, not just approximate. 

4.2. The set-up. Let A be a complex Perron number of degree d. Let 
q[x] G 7j[x] be the minimal polynomial for A; it is monic and of degree 
d. Furthermore the roots Ao,Ai,... , A^-i of q(x) satisfy |Aj| < |A| 
except for the roots A and A. The roots of q(x) are called the Galois 
conjugates of A. Assume A = A. 

Let K = Q[A]. For each % with < % < d — 1 there is an embedding 
Oi : K — > C, sending an element p(X) (a polynomial in A with ratio- 
nal coefficients) to p(Aj). If A, is real the image of this embedding is 
contained in R. 

Let r be the number of real Galois conjugates, and 2c the number of 
non-real Galois conjugates of A, so that r + 2c = d. 

There is a natural embedding a: K — > W x C c = R d which is the 
product of these individual embeddings Oi (taking only one embedding 
for each complex conjugate pair). 

At this point we must distinguish between the case A real and the 
case A is nonreal. 

If A is nonreal, let W = M. d ; let 7r be the linear projection from W to C 
such that Vx G K, ir<r(x) = x. Let m x be the linear map <r(K) — > cr(K) 
induced by multiplication by A in K, i.e. m\{a{x)) = a(Xx). Then m\ 
extends to a linear map of M. d whose eigenvalues are exactly the \, 
and eigenspaces are given by the coordinate axes/planes. Let A = O 
denote the ring of algebraic integers in K. The image <r(A) is a discrete 
lattice in W invariant under m\. 
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If A is real let W = R d x R d and it be the projection tt: W — > M 2 
satisfying 7r(a(x) x = (x, y). Let m-A be the linear map on W 

induced by multiplication by A on K. The eigenvalues of m\ are the 
Aj, each with multiplicity 2. Let i = (!) x where O is the ring of 
algebraic integers in K. The image o~(A) is again discrete lattice in W 
invariant under m\. 

In either case let V\ C W be the eigenplane (2-dimensional eigenspace) 
for mA corresponding to the eigenvalue A. Then the map tt: W — > M 2 
can be thought of as projection along eigenvectors of W onto V\. 

We use the Euclidean metric on both M 2 and W. If x G M 2 or x G 
let Nr{x) denote the neighborhood of radius R around x in M? or VF, 
respectively. Similarly for a set 5 C IR 2 or S C let N R (S) denote 
the ^-neighborhood of S. For v G K let denote the distance from 
cr(v) to the plane V\. This is the length of the "vertical" component of 
a{v). 

For a triangle t C M 2 with vertices fi, t>2, ^3 G -R', let a(t) denote the 
triangle in W with vertices a{v\),a{v2),a{v^). 

Lemma 5. Given a triangle fcl 2 with vertices v\, t>2, t>3 G .4 and an 
e > 0, there is an integer n such that a(\ n t) is almost parallel to V\, 
that is, the projection tt is a (1 + e)-biLipschitz mapping from a(X n t) 
to its image X n t. 

Proof. A non-zero vector v G W lying in the plane of o~(t) has a 
non-zero component in the direction of the plane V\. Since A is the 
strictly largest eigenvalue in modulus of m\, for n sufficiently large the 
component in direction V\ of m™v dominates all the others, and so 
rn\{v) lies almost parallel to V\. The result follows. □ 

4.3. The actors. Let T C M 2 be a polygon with vertices in A, having 
the origin in its interior and sufficiently "round" so that AT contains 
T in its interior. Triangulate the annulus AT — T with triangles whose 
vertices are in A (see Figure H). 

Let M be an integer large enough so that: 

1. M > 2|A| 2 |||f i — III f° r an Y two vertices v i, Vi in the triangulation 
of AT - T Q . 

2. M > |A| 2 |t>i — U2I for any two such vertices. 

3. For any x G W, Nm/(2\\\){%) contains a point of the discrete lattice 
a(A). 

Let T be the set of triangles t = (v±,V2, v 3) with noncollinear vertices 
vx,V2,V3 G A satisfying \\\vi — Vj\\\ < 3M and \vi — Vj\ < 3M for all 
i,j G {1,2,3}. For any t G T, o~(t) is a triangle with vertices in the 
discrete lattice a (A) and edges of bounded length. Thus T is a finite set 



8 



RICHARD KENYON 



^T : 




Figure 2. The triangulation of AT . 

up to translation, i.e. if we consider (vi, v 2 , v 3 ) = (v\ +a,v 2 + a, v 3 + a) 
for all Vi,V2, t>3, a G A. 

Let X be the set of surroundings of elements of T by elements of 
T. By definition a surrounding of t G T is a collection X of elements of 
T, with t G X, which tiles in an edge-to-edge fashion a neighborhood 
of t and such that every element of X touches t in at least a vertex (see 
Figure §.) Note that X is also a finite set up to translation. 




Figure 3. A surrounding. 

Recall that for t = (vi,V2,v 3 ) G T by a{t) we mean the triangle in 
W with vertices (a(vi), <t(w 2 ), a(v 3 )). For X & X, by &(X) we mean 
the union of the <r(t)'s for each triangle t in X. So cr(X) is a piecewise 
flat surface in W . 

Let 9 > be the smallest vertex angle of any triangle in T. Let 
T 2 = 2M and n = r 2 / 1 sin{9 /2). Then in any homothetically scaled 
copy of a triangle in T, the r 2 -neighborhoods of any two edges do not 
intersect outside of the ^-neighborhoods of the vertices (Figure f|). 
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Figure 4. Neighborhood of a vertex. 

Let n be an integer large enough so that for each t G T: 

1. X n t is 1/2- flat in the sense of Lemma [S] (the projection 7r from 
a(X n t) is 3/2-biLipschitz). 

2. the inscribed radius of X n t is at least 2diam(AT ) + 2r 2 . 

4.4. subdividing an element of X. Let I G ^ be a surrounding 
of t G T. Each triangle in X n X is large and flat by our choice of n. 
The subdivision proceeds in three stages: first we must subdivide a 
neighborhood of the vertices of X n t, then a neighborhood of the edges 
of X n t, and then X n t itself. (In the final construction we'll use n + 2 in 
place of n; the rules are the same, though.) 

4.4.1. near a vertex. Let X G X be a surrounding of t G T. Let t> be 
a vertex of t Let iV = iV^A^) C X n X. Let 

F„ = {y G -4 n iV : - A n X||| < M/|A|}. 

Here by \\\y — A n X||| we mean the distance vertically from cr(y) to 
<r(X n X). This is the Euclidean distance d(a(y),z) where z G W is 
defined by ir(z) = y and z G <r(A n X). Then Y v is a discrete subset of 
N . (Later on Y v will be the vertices in a triangulation of N in which 
the triangles are in T.) 

4.4.2. near an edge. Let e = = IJjVj be an edge of t and let 

AT' = iV r2 (A n e) - (N ri (X n Vi) U X ri (A n ^)) 

be the r 2 -neighborhood of X n e which is outside the vertex-neighborhoods 
defined in the previous paragraph. Let 

Y e = {yeAnN': \\\y-X n X\\\<M/\X\}, 

where \\\y — A n X||| is as defined above. Again, Y e is a discrete subset of 
N' . 
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4.4.3. for a triangle. Once the sets Y v ,Y e have been defined as above 
for each vertex and edge of t, define 

Y t = {y G AH X n t : |||y-A n X||| < M/\X\}. 

Let 

Y = Y t U Y Vl U Y V2 U Y V3 U Y ei2 U Y e23 U Y e31 ; 

then Y is a discrete set of points in the r 2 -neighborhood of X n t. 

Triangulate the set Y using the Delauney triangulation, leaving out 
triangles whose circumcircle is not completely contained in N r2 (X n t). 
Each triangle t! = (wi,w 2 ,w 3 ) in this triangulation satisfies: 

\wi - Wj\ < M/\X\ (1) 

since the circumcircle of t' contains no point of Y, and so must have 
diameter bounded by M/\X\ by condition (3) on M. Furthermore each 
triangle satisfies 

1 Ml 
\\wi - Wj\\ < \\\wi - X n X\\\ + \\\X n X - Wj\\\ + -\wi - Wj\ < — (1 + 1 + -) 

2 ' |A| fa 

here using the fact that a(X n X) is |-flat (condition (1) on n). 

Since |A| > 1, these imply t' G T. We have defined a subdivision 
of a neighborhood of t into triangles in T. In addition, for each such 
triangle t', we have At' G T by (0) and (0). 

If two surroundings X±,X 2 (of ti, ^ respectively) overlap, i.e. t 2 G Xl 
and ti G then their subdivisions agree near where X n ti intersects 
X n t 2 (by the local nature of the Delauney triangulation), so that XiUX 2 
has a well-defined subdivision. 



4.4.4. the central tile. We need to alter slightly this subdivision rule 
for triangles so as to insert a copy of the region XTq in the subdivision. 
Let t and X be as before with Y v , Y e defined for each vertex and edge. 
Put a translate C of XTq in the interior of X n t so that for any vertex w 
of C, we have 

HI it; - A n t||| < M, (3) 

and C does not intersect the sets N r2 (X n e) for edges e of t. This can 
be accomplished by conditions (1) on M and (2) on n: just put C near 
the center of the inscribed circle of X n t, so that vertices of C are in A, 
and so that <j(C) is (vertically) within M of a(X n t). Take 

Y' = {y G Y | y £ interior (C)}. 
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As before use the Delauney triangulation of Y', except near C: the 
edges of the boundary of C may unfortunately not be edges of De- 
launey triangles. Take only those Delauney triangles whose circumcir- 
cle doesn't penetrate into C. The remaining untriangulated area is an 
annulus around C, and for any vertex x of the outer boundary Z of 
this annulus, the distance to the inner boundary dC is at most M/\X\, 
since x is on a circle of diameter at most M/|A| which touches C. In 
particular x is at distance at most 2M/|A| from a vertex of dC, since 
the edges of dC have length at most M/\\\ 2 < M/\X\. So we can tri- 
angulate this annulus in the following manner: for each edge C\C2 of 
dC, connect C\,C2 to their closest common neighbor on Z. Connect 
the remaining vertices C3 on Z to the only place now possible on DC 
which doesn't cross edges already present. Each such triangle has edge 
lengths at most 3M/|A|. By the condition (^) on T, each triangle is in 
T. 

This defines a subdivision of a neighborhood of t into triangles in T 
and a region C which is tiled by a translate of the original tiling of XTq. 
Furthermore triangles t near the edges of \ n t satisfy: \t e T. (Figure 



I) 




Figure 5. Subdividing a triangle. 



4.5. the tiling. We can now define the tiling of the plane as follows. 
Start with a tile T at the origin. When we multiply by A, the image 
ATq subdivides into T union the given triangulation of AT — T . We 
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multiply by A n+2 more times, mapping the annuli AT — T successively 
to A 2 T - AT , A 3 T - A 2 T , and so on, up to A" +3 T - A n+2 T . At this 
stage we have a concentric sequence of annuli (each with a triangulation 
homothetic to the triangulation in ATq), and a tile To at the origin. 

Each triangle except those in the innermost and outermost annulus 
has a surrounding by triangles. The surroundings of triangles in the 
second innermost annulus are in X, since the triangles in the first three 
annuli are in T by conditions (1) and (2) on M. So each surrounding is 
a homothetic copy of an element of X. At this point the surroundings 
of the triangles in the 2nd outermost annulus are of the form A n A for 
some X £ X. Now subdivide (as in section |4.4.1|j4T4"Z?D all vertices 
and edges which are on the boundary between the outermost and 2nd 
outermost annulus (i.e. on \ n+2 dT ). 

Multiply once more by A. Subdivide all triangles in the 2nd out- 
ermost annulus A n+3 T — A ra+2 T using the subdivision as in sections 



4.4.1| , f4.4.2| , [4.4.4| . That is, first subdivide any vertices which have not 
yet been subdivided (i.e. were not subdivided at the previous stage), 
then any edges which have not yet been subdivided, then finally subdi- 
vide triangles when all their vertices and edges have been subdivided. 
Those triangles adjacent to X n+2 dT have a previously defined subdivi- 
sion on one or more of their vertices and edges. The new subdivisions 
will by construction agree with the old subdivisions near these vertices 
and edges. Note that in this way the subdivision of a triangle depends 
only on the triangles in its surrounding (i.e. on the appropriate element 
of X). 

We now continue multiplying by A and subdividing, with the fol- 
lowing rule: label the triangles by their "size" : those in XTq — Tq are 
labeled 1, the image of a triangle of label k is a triangle of label k + 1 
unless k = n + 2, in which case the image is subdivided into a copy of 
T and lots of small triangles in T of label 1 again. 

As we continue multiplying these tilings fill out the whole plane, and 
there is a unique limiting tiling since larger and larger portions around 
the origin are fixed. 

Adjacent triangles in the limiting tiling have the same size or else 
their size differs by 1 mod n + 2. If adjacent triangles have size both 
< n + 2, then they don't subdivide on multiplication by A. If the sizes 
are n + 2, n + 2 then their subdivisions agree by definition. If the sizes 
are n + 1, n + 2, then the n + 2-sized tile subdivides, so the common 
egde or vertex subdivides. At the next multiplication by A, the other 
tile will subdivide with a subdivision which by construction matches 
the existing subdivision along the common edge or vertex. If adjacent 
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sizes are 1 and n + 2, then they just came from a pair of tiles of sizes 
n + 2 and n + 1, so again their subdivisions agree. 

4.6. refining the boundaries. We have defined a subdivision rule 
for triangles, or rather, elements of X, and a tiling of the plane with 
scaled elements of X U {AT }. To make this a true self-similar tiling 
we need to redraw the boundaries of the triangles and of T so that the 
image of a tile covers exactly a set of tiles in the tiling. 

Let X G X be a surrounding of t G T. For each edge e in i (between 
vertices Vi and t>j of t) draw a polygonal arc a e from X n+2 Vi to A n+2 t>j, 
which is contained in iV r2 (A n+2 e), and has edges contained in the set 
of edges of triangles in the triangulation of N r , 2 (\ n+2 t). 

The part of the path a e in N ri (X n+2 Vi) should depend only on this 
neighborhood; similarly the part of the path in N ri (X n+2 Vj) should 
depend only on that neighborhood, and the part of the path along the 
edge should depend only on N r2 (X n+2 e). This is to guarantee that the 
paths for different surroundings agree whenever there is overlap. For 
example, if triangles t\ and t 2 are adjacent along an edge e, then the 
path a e along X n+2 e will be the same for the surrounding of t\ as for 
t 2 . 

Furthermore the choice of paths in vertex neighborhoods N ri (X n+2 Vi) 
should be made so that the paths going out along different edges inci- 
dent at the vertex should be non-crossing and diverging (although 
such paths may have common initial segments, they are disjoint except 
for a connected initial segment). See Figure || 




Figure 6. Choosing a path around t. 

Define j3 e = X~ n ~ 2 a e , so that (3 e is a polygonal arc between vertices 
of t which follows closely the edge e. Define 7* to be the concatenation 
of the 3 edge paths /3 ei2 , /3 e23 , /3 e31 , removing any "backtracking" which 
may occur near the three vertices: that is, if the last edge of (3 ev2 is 
the same as the first edge of /3 e23 in the reverse direction, remove these 
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two edges from 7^, and so on. So 7/ is a simple closed polygonal path 
running once around t. 

The paths \ n+2 ^l are composed of edges of triangles in T. To each 
such edge ej is associated its arc (3 ej . Define 7 2 by replacing each edge 
\~ n ~ 2 ej in 7^ with a copy of X~ n ~ 2 j3 ej , again removing backtracks near 
the vertices of the ej. Inductively define 7 t fc+1 by replacing each edge 
\~( n+2 > k ej in 7 4 fc by a copy of \~( n+2 ' k /3 e . and removing backtracks. 

These paths 7^ converge as k — >• 00 to a simple closed curve 74 
running around t: 7^ lies in the cA~ fe ( n+2 ) -neighborhood of 7 4 fc_1 , for 
some constant c, so the paths converge in the Hausdorff metric. The 
limits are arcs by our choice of n, which is so large that any poten- 
tial non-injectivity in 7^ must occur near a vertex of 7f fc_1 , where by 
construction it cannot occur. 

For each X E X, the Jordan curve 7* forms the boundary of a tile 
Tx- A similar construction gives a Jordan curve for the boundary of 
T . These tiles T x have the property that \ n+2 T x subdivides exactly 
into other tiles in {T } U {T x > '■ X G X}. Also AT subdivides exactly 
into a copy of T and certain T x . 



4.7. proof of self-similarity. Let S = {T } U S x U £ 2 U . . . U S n+2 , 
where Si = {\ l T x \ X G X}. We have now constructed a tiling of 
the plane, with tiles which are translates of tiles in a subset of S (not 
every element of S may actually occur in the tiling). By construction 
when we multiply the tiling by A and subdivide each tile according to 
the rules established above (T subdivides into the tiles as defined in 
section |2|, each tile in St maps homothetically to a tile in Sj+i unless 
i = n + 2 in which tile subdivides as described in section |4.4.4|) . 

Furthermore tiles of the same type (i.e. those arising from the same 
element of X, and which are of the same size) subdivide in the same 
way. 

Furthermore each tile occurring in the tiling arises from eventual 
subdivision of the tile T , since T has the origin in its interior. Each 
tile eventually has a copy of T in its subdivision. So the subdivision 
matrix is primitive, and hence the tiling is quasiperiodic (Lemma|). 

This completes the proof. 



5. Tilings invariant under rotations 

Let T be a self-similar tiling invariant under e %e and having expansion 
constant A. If A is real, then by redefining the subdivision rule we have 
that T is self-similar with expansion e %6 \. Thus we can assume A is not 
real unless e %e = — 1 . 
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In IJTTJ it is shown that for any self-similar tiling with expansion A ^ R 
there is a homothetic copy K' of Q[A] such that translations between 
tiles of the same type must be in K' . So if the tiling is invariant under 
x — > xe ld , then if a ^ is a translation between tiles of the same type, 
both a and ae ld are in K', so their ratio e ld must be in 

In case e %e = — 1 we also have e %e ( 



To construct a tiling invariant under rotation by 9 = 27r/m, modify 
the construction of section |] as follows. 

Select the tile To to be invariant under e ld . Then the region between 
XTq and To has a e l£l -invariant triangulation. 

The construction now proceeds as before. Note that e ld is an al- 
gebraic integer, and multiplication by e ld acts as an isometry on the 
space W and the lattice <j(A). Since the subdivision defined in section 
Q| depends solely on the metric properties of the objects involved, 



the subdivision is natural in the sense that two elements X 1 ,X 2 E X 
(surroundings of t\ , t 2 respectively) which are isometric copies of each 
other have subdivisions which are isometric. (Note that even if X\ 
and X 2 overlap, the naturality of the subdivision implies that they will 
subdivide isometrically) . 

The only place in the construction where we again had a choice was 
in defining the paths running around the tiles. However we claim that 
we can make this choice in a way which only depends on the isometry 
type of the relevant vertex or edge. The only difficulty arises when 
isometric surroundings overlap. If ti and t 2 meet at a vertex v and 
the isometry from t\ to t 2 takes v to v, then the subdivision near v 
is already ^-invariant, so one can choose paths running out from that 
vertex in a ^-invariant fashion. Similarly if t± and t 2 meet along an edge 
e, and the isometry taking t\ to t 2 takes e to e, then the subdivision of 
e is invariant under this isometry and so one can choose the path a e 
to also be invariant. These choices imply the naturality of the entire 
construction. 

So the resulting tiling will be invariant under e %e . 



6. Construction from free group endomorphisms 

In this section we give a construction for self-similar tilings whose 
expansion coefficient is a complex Perron number which satisfies 

A" - pX 1 - 1 + q\ + r = (4) 

for some n > 3 and integers p, q > 0, r > 1. 

We will first do the case n = 3. Let a, b, c be vectors pointing in 
different directions in R 2 . Let F be the set of compact polygonal paths 
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starting at the origin, each of whose edges is a translate of ±a, ±6, or 
±c, and which are non-backtracking in the sense that along the path a 
segment x is never immediately followed by —x. 

Thus to each element in F is associated a unique element in the 
free group on three symbols F(a, b, c). Indeed, one can easily define a 
product on F in such a way that the map F(a, b,c) —>■ F is a group 
isomorphism: to obtain the product of two paths, translate the second 
to the end of the first, and then cancel any "backtrackings". 

Let /: F(a, b,c) — > F be this isomorphism. 

Define an endomorphism 0: F(a, b, c) — > F(a, b, c) by: 

0(a) = b 

m = c , 

0( c ) = d>a~ r b- q . 

Consider the three commutators [a, b] = a6a _1 6 _1 , [b, c], and [a, c]; 
they represent three closed paths. Their images under can be written: 

0([a,o]) = [0(a),0(6)] = M (5) 
<f>([b,c]) = [c,c p a- r b^] =c p a- r [a r ,c]{b- q [b\c]b q )a r c- p (6) 
0([a,c]) = [&,c?a- r &- 9 ] = [6,c p ]c?a- r [a r ,6]a r c- p (7) 

Using the identities 

[x n ,y] = {x n - l [x,y]x 1 - n ) ■ ■ ■ {x[x,y]x- l )[x,y} 

and 

[x,y n ] = [x,y](y[x,y]y- 1 )...(y n - 1 [x,y]y 1 - n ) 

in the above, we see that each of 0([a, b]), <f)([b, c]) and 0([a, c]) can be 

written as a product of conjugates of [a, ft], [6, c], and [a, c]. Interpreting 

this in terms of paths in M 2 , the closed paths f<f)([a, b]), f<fi([b, c]), f<f>([a, c]) 

can each be tiled by translates of copies of the parallelograms f([a,b)),f([b,c]),f([a 

Geometrically, the images f([a,b]),f([b,c)),f([c,a]) are "archtiles" 
for a self-similar tiling of the plane. Let S = {[a, b], [b,c], [a,c]}. The 
images f(<p(x)), where x G S can be tiled by translates of tiles in S: 
the exact translates can be calculated from the formulas in 

Take a, b, c to be the vectors 1, A, and A 2 6 C respectively. The 
equation (|]) implies that the argument of A is less than tt/2 degrees, 
so that f[a,b], f[b,c], f[a,c] each have the same orientation. (In par- 
ticular this implies that the tiling will be non-overlapping.) Then 
/0(a) = /(&) = Xf(a), f<f>(b) = /(c) = A/(6), and /0(c) is a path 
whose endpoint is at A/(c). 
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Now the closed paths f(<j) n (x)) for ieS converge after rescalingj] to 
the boundaries of certain tiles Ti, T 2 , T3, and ATj can be tiled by exactly 
by translates of Ti,T 2 and T3. This is sufficient to make a self-similar 
tiling of the plane. 

An example is given in Figures [7|,|8|J| and |I[ 




Figure 7. The case (p,q,r) = (1,2,1) (A « .696 + 
1.436z). The three archtiles /([a, 6]), /([6, c]), /([a, c]). 




Figure 8. The successive approximations to the subdi- 
vision of [b, c]: here we took c]) for n = 2, 3, 4, 5, 6, 
and 10. 




Figure 9. The subdivisions of tiles Ti,T 2j T 3 : 7\ gives 
a T 2 , T 2 gives two copies of T 2 and a copy of T 3 , and T 3 
subdivides into a T 2 and a Ti. 

For the case n > 3, a similar construction will work. Here we 
must use the free group on n letters with an endomorphism 

l a. rigorous argument for convergence can be made along the lines of Dekking 

m 
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4>(a,i) = a i+ i for 1 < i < n — 1 and (f>(a n ) = a^a^a^ '. The archtiles 
are the (n 2 — n)/2 words/paths [a*, %] with i < j. 

This same construction in fact works with any endomorphism qb 
which has the property that each 0([aj, a,-]) can be written as a product 
of conjugates of the [a^a,/]. Since conjugates of the [aj,Oj] generate 
the commutator subgroup [F, F] , this condition amounts to the require- 
ment that the conjugates which appear have nonnegative exponent. 

For an arbitrary endomorphism, one can also attempt to change the 
basis of F so that it has the correct form; a necessary condition is that 
the induced linear map on [F, F]/[F,[F, F]] = Z^" n ^ 2 have matrix 
with nonnegative coefficients. 
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